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Abstract 

Methods for the computation of invariants and symmetries of nonlinear evolution, 
wave, and lattice equations are presented. The algorithms are based on dimensional 
analysis, and can be implemented in any symbolic language, such as Mathematica. 
Invariants and symmetries are shown for several well-known equations. 

Our Mathematica package allows one to automatically compute invariants and sym- 
metries. Applied to systems with parameters, the package determines the conditions on 
these parameters so that a sequence of invariants or symmetries exists. The software 
can thus be used to test the integrability of model equations for wave phenomena. 

1 The Key Concept: Scaling Invariance 

The ubiquitous Korteweg-de Vries (KdV) equation from soliton theory, 

(1) ut = 6uu x + u 3x , 

is invariant under the dilation (scaling) symmetry (t,x,u) — > (A _3 i, A _1 x, A 2 u), where A is 
an arbitrary parameter. Obviously, u corresponds to two derivatives in x, i.e. u ~ d 2 /dx 2 . 
Introducing weights, w(u) = 2 if we set w(d/dx) = 1. Similarly, d/dt ~ 9 3 /cAr 3 , thus 
w{d/dt) = 3. The rank of a monomial equals the sum of all of its weights. Observe that 
(H) is uniform in rank since all the terms have rank R = 5. 

Likewise, the Volterra lattice, which is one of the discretizations of ([j), 

(2) il n = U n (U n+ l - U n _l), 

is invariant under (i, u n ) — > (A -1 i, Xu n ). So, u n ~ d/dt, or w(u n ) = 1 if we set w(d/dt) = 1. 
Every term in (|2|) has rank R = 2, thus @ is uniform in rank. 

Scaling invariance, which is a special Lie-point symmetry, is common to many integrable 
nonlinear partial differential equations (PDEs) such as (|l]), and nonlinear differential- 
difference equations (DDEs) like (|2|). Both equations have infinitely many polynomial 
invariants pi ^] and symmetries Q. In this paper we show how to use the scaling invariance 
to explicitly compute polynomial invariants and symmetries of PDEs and DDEs. 



2 Computation of Invariants 

For PDE systems as in Table 1, the conservation law Dtp = D X J connects the invariant 
(conserved density) p and the associated flux J. As usual, Dj and D x are total derivatives. 
Most polynomial density- flux pairs only depend on u, u x , etc. (not explicitly on t and x). 
Integration of the conservation law with respect to x yields that P = p dx is constant 
in time, provided J vanishes at infinity. P is a conserved quantity. 
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The first three (of infinitely many) conservation laws for ([!]) are 

(3) D t («) = D,, (3u 2 + u 2x ), D t (u 2 ) = B x (4u 3 - u\ + 2^ 2:E ) , 
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(4) Dt(u 3 - -u 2 ,) = B x (-u 4 - 6uu 2 x + 3u 2 u 2x + ~ u x u 3x ). 

The densities p = u,u 2 ,u 3 — \u\ have ranks 2,4 and 6, respectively. 
Conserved densities of PDEs like (|l]) can be computed as follows: 

• Require that each equation in the PDE system is uniform in rank. Solve the resulting 
linear system to determine the weights of the dependent variables. For (|l|), solve 
w(u) + w(d/dt) = 2w(u) + 1 = w(u) + 3, to get w(u) = 2 and w(d/dt) = 3. 

• Select the rank R of p, say, R = 6. Make a linear combination of all the monomials 
in the components of u and their x-derivatives that have rank R. Remove 'equivalent' 
monomials, that is, those that are total ^-derivatives (like u^ x ) or differ by a total x- 
derivative. For example, uu 2x and u 2 x are equivalent since uu 2x = \ {u 2 ) 2x — u 2 x . For (|l]), 
one gets p = ciu 3 + c 2 u 2 x of rank R = 6. 

• Substitute p into the conservation law, eliminate all ^derivatives, and require that the 
resulting expression is a total x-derivative. Apply the Euler operator [|J to avoid integration 
by parts. The remaining part must vanish identically. This yields a linear system for the 
constants q. Solve the system. For (|l|), one gets c\ = l,c 2 = —1/2. 

See Q for the complete algorithm and its implementation. See [|| for an integrated 
Mathematica Package that computes invariants (and also symmetries) of PDEs and DDEs. 



Table 1 
Invariants and Symmetries 





Continuous Case (PDEs) 


Semi-discrete Case (DDEs) 


System 


u t = F(u, u x ,u 2x , ...) 


U n = F(..., U„_l, U n , U n+ i, ...) 


Cons. Law 


B t p = D X J 


Pn — Jn Jn+1 


Symmetry 


D t G = F'(u)[G] = fF(u + eG)| e=0 


D t G = F'(u n )[G] = fF(u n + eG)| e=0 



For DDEs like (|2|), compute the weights in a similar way. Determine all monomials of 
rank R in the components of u n and their ^derivatives. Use the DDE to replace all the 
i-derivatives. Monomials are 'equivalent' if they belong to the same equivalence class of 
shifted monomials. Keep only the main representatives (centered at n) of the various 
classes. Combine these representatives linearly with coefficients q, and substitute the form 
of p n into the conservation law p n = J n — J n +±. Remove all t-derivatives and pattern-match 
the resulting expression with J n — J n +±. Set the non- matching part equal to zero, and solve 
the linear system for the c\. Determine J n from the pattern J n — J ra +i. For (||), the first 
three (of infinitely many) densities p n are listed in Table 2. Details about the algorithm 
and its implementation are in ||, ^[ |3| . 



3 Computation of Symmetries 

As summarized in Table 1, G(x, t, u, u x , U2 X , •••) is a symmetry of a PDE system iff it leaves 
it invariant for the change u — > u + eG within order e. Hence, D^(u + eG) = F(u + eG) 
must hold up to order e. Thus, G must satisfy the linearized equation D t G = F'(u)[G], 
where F' is the Frechet derivative: F'(u)[G] = t|F(u + eG)| £=0 . 

Using the dilation invariance, generalized symmetries G can be computed as follows: 

• Determine the weights of the dependent variables as in Section 2. 

• Select the rank R of the symmetry. Make a linear combination of all the monomials 
involving u and its x-derivatives of rank R. For example, for (||), G = c\ u 2 u x + 02 u x U2 X + 
C3 uv,3 X + C4 u§ x is the generalized symmetry of rank R = 7. 

• Compute D^G. Use the PDE system to remove all i-derivatives. Equate the result to the 
Frechet derivative F'(u)[G]. Treat the different monomial terms in u and its x-derivatives 
as independent, to get the linear system for q. Solve that system. For (|l|), one obtains 

(5) G = 30u 2 u x + 20u x U2 X + Wuu 3x + u^ x . 

The symmetries of rank 3, 5, and 7 are listed in Table 2. They are the first 
three of infinitely many. For DDEs like (|2|), G(..., u n _i, u n , u n+ i, ...) is a symmetry iff 

Table 2 
Prototypical Examples 





Korteweg-de Vries Equation 


Volterra Lattice 


Equation 


u t = 6uu x + u 3x 


U n = U n (u n+ i - Un-i) 


Invariants 


p = u p = u 2 
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P = U ~ 2 U x 


Pn = Un Pn = U n (\u n + U n+ i) 
Pn = \Un + U n U n+ i (u n + U n+ i + U n+2 ) 


Symmetries 


G = u x G = 6uu x + us x 
G= 30u 2 u x + 20u x U2 X + Wuu3 X +u§ x 


G = u n u n+ i(u n + U n+ l + U n+ 2) 
-U n -iU n (u n -2 + n n _i + u n ) 



the infinitesimal transformation u ra — > u n + eG leaves the DDE invariant within order 
e. Consequently, G must satisfy 4^- = F'(u„)[G], where F' is the Frechet derivative, 

F'K)[G] = |F(u n + eG)| e=0 . 

Algorithmically, one performs the following steps: First compute the weights of the 
variables in the DDE. Determine all monomials of rank R in the components of u n and their 
t-derivatives. Use the DDE to replace all the i-derivatives. Make a linear combination of 
the resulting monomials with coefficients Cj. Compute DtG and remove all u n _i, u n , Un+i) 
etc. Equate the resulting expression to the Frechet derivative F'(u n )[G] and solve the 
system for the Cj, treating the monomials in u n and its shifts as independent. Details are 
in |Q, For (0), the symmetry G of rank R = 3 is listed in Table 2. There are infinitely 
many symmetries, all with different ranks. 

See for the complete algorithm and its implementation in Mathematica, and j|] for 
an integrated Mathematica Package that computes symmetries of PDEs and DDEs. 



Notes: 

(i) If PDEs or DDEs are of second or higher order in t, like the Boussinesq equation in [|J , 
we assume that they can be recast in the form given in Table 1. 

(ii) A slight modification of the methods in Section 2 and 3 allows one to find invariants 
and symmetries that explicitly depend on t and x. See next section for an example. 

(iii) Applied to systems with free parameters, the linear system for the c, will depend on 
these parameters. A careful analysis of the eliminant leads to conditions on these parameters 
so that a sequence of invariants or symmetries exists. 

(iv) For equations that lack uniformity in rank, we introduce one or more auxiliary 
(constant) parameters with weights. After the form of the invariant or symmetry is 
determined, the auxiliary parameters can be reset to one. 

(v) Higher-order symmetries, such as (0) lead to new integrable evolution equations. For 
example, Ut = 30u 2 u x + 2§u x u% x + Wuus x + u^ x is the completely integrable fifth-order 
KdV equation due to Lax. 

Details about these 5 items are given in (l|, ||, ||, |j| . 
4 Examples 

4.1 Vector Modified KdV Equation 

In 0, Eq. (4)], Verheest investigates the integrability of a vector form of the modified KdV 
equation (vmKdV), which upon projection, reads 

ut + 3u u x + v u x + 2uvv x + u 3x = 0, 

(6) v t + 3v 2 v x + u 2 v x + 2uvu x + v 3x = 0. 

With our software InvariantsSymmetries.m [|j we computed the following invariants: 

(7) pi = u, p 2 = v, P3 = u 2 + v 2 , 

(8) p 4 = l -(u 2 + v 2 ) 2 -(u 2 x + vl), 

(9) p 5 = ^x(u 2 + v 2 )-^t(u 2 + v 2 ) 2 + t(u 2 x + v 2 x ). 

Note that the latter invariant depends explicitly on x and t. Verheest 0j has shown that 
(0) is non-integrable for it lacks a bi-Hamiltonian structure and recursion operator. We 
were unable to find additional polynomial invariants. Polynomial higher-order symmetries 
for (||) do not appear to exist. 

4.2 Extended Lotka-Volterra Equation 

Itoh U studied the following extended version of the Lotka-Volterra equation (||) , 

fc-1 

(10) u n = Y^(u n _ r - u n+r )u n . 

r=l 

For k = 2, ( |l0|) is @, for which three invariants and one symmetry are listed in Table 2. 

In Q, we give two additional invariants; in |g] we list two more symmetries. 

For (|io|), we computed 5 invariants and 2 higher-order symmetries for k = 3 through k = 5. 



Here is a partial list of our results: 

Case 1: k = 3 

Invariants: 

(11) Pi = Un, P2 = ^Un + U n (u n+ i+U n+2 ), 

P3 = + U 2 n (u n+ i + U n+2 ) + U n {u n+ i + U n+2 ) 2 

(12) +u n (u n+ iu n+3 + U n+2 U n+3 + u n+2Un+4 J • 

Higher-order symmetry: 

G = u n (u n +i + u n+2 - u n - 2 - u n -i) + u n [(u n+ i + u n+2 ) 2 - (u n _ 2 + u n -i) 2 ] 

(13) +u n [u n+ iu n+3 + u n+2 u n+3 + u n+2 u n+4 - (-u n _4-u n _2 + ■u n _ 3 n„_ 2 + U„_ 3 U„_i)]. 

Case 2: = 4 

Invariants: 

(14) pi = v^, p 2 = ]^u 2 n + u n {u n+ i + u n+2 + u n+ s), 

P3 = + U n (u n+ 1 + U n+2 + U n+3 ) + U n (u n+ \ + U n+2 + U n+3 ) 2 

(15) +n n (n n+ i-u n+4 + -u n+2 u n+4 + u n+3 u n +4 + u n+2 u n+5 + u n+3 u n+5 + u n+3 u n+6 ). 
Higher-order symmetry: 

G = U n [u n+ lU n+ 4 + U n+2 U n+ 4 + U n+3 U n+ 4 + U n+2 U n+ 5 + U n+3 U n+ ^ + U n+3 U n+ Q 
-{Un-6Un~3 + W„_ 5 U n _ 3 + U n _ 4 U n „ 3 + U n _ 5 M n _ 2 - U n -4U n _ 2 + M n _ 4 U n _i)] 
+M|j[K+l + Un+2 + Un+3) 2 ~ U n (u n -3 + U n - 2 + U n _i) 2 ] 

(16) +u 2 [u n+ i + U n+2 + U n+3 - (u n - 3 + n n _ 2 + U n -i)]. 
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